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In this paper we study the footprint of cosmic string as the topological defects in the 
very early universe on the cosmic microwave background radiation. We develop the method 
of level crossing analysis in the context of the well-known Kaiser-Stebbins phenomenon for 
exploring the signature of cosmic strings. We simulate a Gaussian map by using the best 
fit parameter given by WMAP-7 and then superimpose cosmic strings effects on it as an 
incoherent and active fluctuations. In order to investigate the capability of our method to 
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detect the cosmic strings for the various values of tension. G/i, a simulated pure Gaussian map 
is compared with that of including cosmic strings. Based on the level crossing analysis, the 
superimposed cosmic string with G/i > 4 x 10 9 in the simulated map without instrumental 

^ ! 

noise and the resolution R = V could be detected. In the presence of anticipated instrumental 
V£2 ■ noise the lower bound increases just up to G/i > 5.8 x 10~ 9 . 
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I. INTRODUCTION 

One of the most interesting predictions of quantum field theory in the domain of cosmology is 
the possibility of transition between different vacuum states during the expansion of the universe. 
Depending upon the topology of these states, a series of stable topological defects such as domain 
walls, monopoles and cosmic strings can be formed. Regarding the formation of strings, they 
may have self-interactions causing the formation of closed loops according to the so-called inter- 



commutation 



10( | . In principle, string loops start to oscillate and emit stochastic gravitational 



waves which results in their annihilation, 



while the infinite straight strings survive up to now. 



Although, in some models of inflation such as false vacuum dominated inflation (in the context of 
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hybrid inflation theories), topological defects are formed at the end of inflation era [12j, but the 
true focus on the issue of cosmic strings was in the recent years, mainly because new models of 
inflation derived from the superstring theory, result in acceptable possibility for the production 



of cosmic strings |!3l4l6l|. This is based on the consideration of cosmic string networks consisting 
of infinite strings, loops, and junctions of two or more strings which will definitely have effects 
such as lensing, CMB polarization, and of course CMB anisotropies. Astrophysical evidence of 
cosmic string strongly depends on the two following parameters, namely: dimensionless string 
tension G\ijc 2 ~ A 2 /M| lanck (G is the Newton's constant and A stands for energy scale when the 
strings are created) and inter-commuting probability, P. The quantity [i shows the mass per unit 
length of cosmic string. The energy density of cosmic string is associated with the scale of phase 
transition and symmetry breaking in which cosmic strings are produced. Generally the formation 
of cosmic strings could occur at very extended ranges of energy scale, e.g. the Grand Unified 
Theory (GUT) scale with A ~ 10 16 GeV corresponding to G/i/c 2 ~ 10 -6 , consequently there is a 
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wide range of Gfi/c 2 s has been supposed |9|, ll7j. It must be noted that the footprints of strings, 
such as anisotropies in CMB are directly affected by their tension. So determining the bounds for 
the tension, directly means limiting the fundamental theory on the basis of which cosmic strings 
are produced. To this end, observing the cosmic string footprints by using various approaches can 
be interpreted as a kind of observational evidence for the low energy limit of the superstring theory 
and would provide the most direct test of string theory and could rule out or constrain on particle 
physics models. 

From theoretical and observational perspectives, there are a dozen constrains on the cosmic 
string's parameters. Recent analysis by using pulsar timing and photometry based on grav itational 
microlensing put a constraint for tension in the range 10 -15 < G[ijc 2 < 10 



COSMOS survey, there is Gfi/c 2 < 3 x 10 7 for cosmic string 



signature of cosmic string wakes has also been explained in Ref. 



reported in 



1MM. Using 



27|]. The 21 cm 



2a ] . Another robust constraint on 



the cosmic string's free parameters comes from the temperature fluctuations at the last scattering 
surface. More recent full analysis and prediction for incoming satellite based surveys can be found in 
Q,y,0, 3l3] • Temperature fluctuations contain the accumulation of anisotropies induced by cosmic 
strings and can be divided into two categories: 1) anisotropies created by the so-called Kaiser- 
Stebbins effect and those related to pre-recombination processes and 2) the stochastic background 



of gravitational waves produced by decaying of string loops 
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321 ] . Calculation of temperature 



angular power spectrum puts the upper bounds G///c 2 < 2 x 10 7 [3ll438j] and G^l/c 2 < 6.4 x 
10~ 7 at 2a confidence interval for Abelian-Higgs case [22|. Nambu-Goto numerical simulation for 
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311 ] . For detecting the anisotropics 



network of cosmic string has been used by A. A. Fraisse et al. 
induced by mentioned strings, temperature gradient magnitude operator has been applied in this 
paper. Ligo and Virgo collaborations recently have reported the newly lower and upper bounds as 
7 x 10 -9 < Gfi/c 2 < 1.5 x 10~ 7 on the stochastic gravitational waves produced by cosmic strings 



39(. The skewness in CMB fluctuations and its dependency on inter-commuting probability has 



been investigated in [40(. Based on B-mode polarization of cosmic microwave background new 



constraint on Gfi has been reported in Ref. 



41|. 



Direct implication based on the explicit recognition of discontinuity in the fluctuation of cosmic 
microwave back gro und radiation is a unique signature of straight cosmic string , namely the Kaiser- 



Stebbins effect 



321 ] . In this part, there are also some upper bounds as well as lower bounds for 



distinguishing cosmic strings. Another method is the so-called Wavelet domain Bayesian denoising 
which has recently been used by D. K. Hammond et al. [421 ] . Since the cosmic string is charac- 
terized by a localized edge like discontinuity, they used a Steerable Pyramid wavelet transform to 
discriminate straight edge from other features. Transformation has been done with 6 orientations 
and 4 spatial scales. According to their results, the signature of cosmic strings in the CMB map 
without noise can be identified for tensions in the range G/i/c 2 > 6. 3 x 10~ 10 while for a noisy 



map the lower bound increases more than one order of magnitude 
reader can also refer to Refs. 
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421 ] . For completing this topic, 



441 ] . The detectability threshold for cosmic string signal in the 



jresence of anticipated noise such as extra-Galactic and Galactic foregrounds has been studied in 



45J. A series of papers published by Brandenberger and his collaborators have been devoted to 



the ability of so-called Cann y a lgorithm to detect the cosmic string and super strings according 
to Kaiser-Stebbins effect 4614481] . They have concluded the lower bound of Gfi/c 2 > 5.5 x 10 -8 . 
Actually in the Canny method, the sharp and strength edges are picked up by means of the mag- 



nitude of the contrast from one side of produced edge to the other one 



461 ] . Non-gaussianity due 



to cosmic strings in addition to other events such as kinetic Sunyaev-Zel'dovich effect have been 
examined in by multi-scale methods (also see [s(J for more discussions). They have argued 
that to find a significant results not only one should use robust methods but also a combination of 
powerful methods is necessary. 

Another most interesting method which is potentially able to discriminate map with and without 
cosmic strings as well as other topological defects is the so-called peak-peak correlation function [51I - 
Qj. In this method, two-point correlation function of the local maxima (or minima) of temperature 
fluctuations on the CMB map at a typical peak threshold is determined. A. F. Heavens et al. have 
demonstrated that this correlation function for a simulated CMB map accompanying cosmic string 
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FIG. 1: A cross-section of space-time perpendicular to a straight cosmic string. Due to deficit angle, the 
light ray passes through the left part of plot to be blueshiftcd with respect to that of received from the right 
part. 



grows up significantly at around 10 — 15 arcminutes, relative to the pure Gaussian map. As 
they claimed, this procedure has no adjustable parameter, consequently it can be powerful for 
test of Gaussianity hypothesis as well as tracking the footprint of other mechanism for producing 
initial fluctuations in the upcoming high resolution observations. Although, they didn't show the 
efficiency of this method for various values of G[i which is the main purpose of this paper, but it 
is interesting reapply the same idea based on hotspots (or coldspots) framework. In our method 
not only we investigate all features of fluctuations instead of local extrema used in the peak-peak 
correlation function, but also the q— moments of relevant quantity are used (see section II for more 
details). 

In this study we are relying on a newly demonstrated method to put a lower limit on the tension 
of straight cosmic strings in the high resolution of CMB map for which we can discriminate their 
signatures. As the effect of cosmic strings, we concentrate on the discontinuities and fluctuations 
in the CMB map arising from the Kaiser-Stebbins effect which is a part of integrated Sachs- Wolfe 
effect and can produce observational consequences on the anisotropies in the CMB map. This 
eature is based on the gravitational lensing which was a well-established phenomena by the time 



32, 



5514581]. Consider a string extended in the direction perpendicular to the paper (see Figure ([I])), 



and two photons emitted from a source in the background. Choose two paths on opposite sides 
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of the string that reach the same observer. If the string is moving in the direction perpendicular 
both to its length and to the line of sight between the source and observer, with the velocity v a , 
a frequency shift will occur between two photons. Thus if we move across the cosmic string when 
sweeping the CMB we can see a jump in the anisotropy map. The amount of this jump is 

ST 

— = 87rG/x|n.(7 s v s x e s )| (1) 

in which n is the direction of observation, v s is the velocity of the string, e s its orientation, and 
7 S is the Lorentz factor for the string. Eq.([T|) can provide a basis for detecting cosmic strings via 
direct observation of the CMB anisotropy map. However, it must be noted that since the Kaiser- 
Stebbins effect may be considerably smaller than the observed temperature fluctuations in WMAP 
data, we cannot expect to detect strings in the WMAP data. Meanwhile we are able to produce 
simulated maps of CMB with any set of parameters. So in order to examine the bounds on the 
string tension imposed from this analysis, we can use simulated maps to apply certain methods of 
statistical analysis on the inhomogeneities to detect the strings or to assume their existence and 
then obtain better bounds on string tension. 

The rest of paper is organized as follows: In Sec. II we introduce level crossing analysis to 
investigate the fluctuation of temperature on the last scattering surface with and without cosmic 
strings. Simulation maps of CMB for various values of so-called string tension using most recent 
observation based on WMAP-7 mission will be given in details in Sec. III. Data analysis by using a 
robust method in complex system from a statistical point of view to distinguish the Gaussian map 
with an without cosmic strings in the presence of expected instrumental noise will be explored in 
Sec. IV. Summary and conclusions are presented in Sec. V. 



II. LEVEL CROSSING ANALYSIS ON THE FLAT SKY MAP 



After the innovation of Rice for calculating the upcrossings and downcrossings of a stochastic 
signal at an arbitrary level 



59l | . many studies have been devoted to signal processing as well as 
fluctuations of height of rough surfaces from the level crossing point of view as a powerful method 
in complex system 



60-65|]. In the level crossing analysis, we are interested in determining the 
number crossings of the temperature fluctuations, N a , at an arbitrary level a. We consider a 
sample function of temperature fluctuations which is defined on the homogeneous and isotropic 
random surface represented by T(n) at direction h relative to the observer placed at the center 
of sphere. For a flat patch of CMB map, we assign temperature fluctuations to each point by 
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FIG. 2: Positive slope crossing at the level T = a. 

T(xi,yi), where X{ and yi demonstrate the coordinate positions. For convenience, the origin of the 
coordinate system of the patch is placed in the left bottom of corresponding map. 

Since we consider the statistical isotropy and homogeneity of temperature random field on the 
surface, without losing generality, we assume a one-dimensional slice of temperature fluctuations 
on a patch of CMB with size G parallel to x or y axes. The size of this signal depends on the 
resolution of data set which is simulated or observed depicted by R = N/Q. Also the numbers of 
pixels will be N 2 . Now consider a sample of an ensemble of one dimensional signal of temperature 
fluctuation, T(k), for which k runs from 1 to N. Suppose n+ denotes the numbers of positive slope 
crossings (upcrossings) of, T(k) — (T) = a, for a typical sample size Q = R x N (see Figure (J2J) ). 
The ensemble averaging for level crossing with positive slope is also given by: 



(2) 



For a statistical isotropic and homogenous fluctuation, the ensemble averaging can be done on 
various one dimensional slices of fluctuations. For convenience, we chose these slices to be parallel 
to x or y axes. Furthermore, the average number of crossings is proportional to the space (time) 



interval G 



60]. Hence: 



iV+(G) oc G 



iV+(G) 



(3) 
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p(y,y') 




FIG. 3: Sketch of joint probability density function of a typical fluctuation and its derivative with respect 
to corresponding dynamical parameter (angle) in the level crossing theory. The shaded area indicates the 
total probability of crossing with positive slope at level y = a. 



which i/+ is the average frequency of positive slope crossing of the level T(k) — (T) = a. It must be 
pointed out that in the long run processes the conservation law for upcrossings and downcrossings 
will be satisfied 60]. In our analysis due to the statistical isotropy and 200 numbers of run over 



ensembles or even more, this conservation law statistically holds and so, we use the upcrossing 
events throughout our analysis. From mathematical and statistical points of view, the frequency 
parameter i/+ can be deduced from the probability distributions associated with the temperature 
fluctuations, T{k) — (T). In order to calculate this probability distribution function, we assume 
two necessary and sufficient conditions for crossing the level T = a by the underlying signal with 
positive slope as follows: 

1) At the beginning of the interval we should have T(fci) — (T) < a. 

2) The slope of signals should be larger or at least equal to the slope of line which is drawn between 
the point at the beginning of interval and the point located on the horizontal line of level a (Figure 
(j2])), namely: 

d[T(k)-(T)] a - [T(h) - (T)] 

de de { ' 



60l465l|. we statistically expect to 



where dO = R. If the above-mentioned conditions are satisfied 
get a high probability for a crossing in the interval dQ. In order to determine whether the above 
conditions are satisfied at any arbitrary location k, we should find how the values of y = T{k) — (T) 
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and y' = are distributed by considering their joint probability density function, p(y,y'). For a 
specific level y = a and interval d@, we focus on the values of y < a and values of y' = (^g) > 
corresponding to the region between the lines y = a and y' = in the plane (y, y') (see Figure 
©)• Subsequently, the probability of positive slope crossing of signal at level y = a in the ciO 
interval is given by: 

/'OO m 

Probability = / dy' I p(y, y')dy (5) 

JO Ja-y'dB 

When d@ — > 0, it is legitimate to put: 

p(y,y') = p{y = a,y') (6) 

Since at large values of y and y', the probability density function approaches zero fast enough, 
Eq.([5|) may be written as: 

roo r-OL 

Probability = / dy \ p(y = a, y')dy (7) 

JO Ja-y'dB 

In this form, the integrand does not depend on y so the first integral can be integrated easily: 

/•oo 

Probability = dO / p(y = a, y')y'dy' (8) 
Jo 

Figure ([3]) shows the area of interest for probability of upcrossings at level a. As mentioned before, 
the average number of crossing with positive slope in interval is ^6. So the average number 
of positive crossings of y = a in interval dQ is equal to the probability of positive crossings of 
y = a in dQ. We can write the average number of upcrossings at the level a in terms of the joint 
probability density function as follows: 



oo 



v£dQ = dQ p(a,y')ydy (9) 



o 



or 



oc 







p(a,y)ydy (10) 
Another useful parameter based on u£ can be introduced as: 

r+oo 

KM= vtW-a\ q da. (11) 



Obviously, for q = the above quantity determines total number of upcrossings for temperature 
fluctuations with positive slope at all levels. For a typical rough fluctuation, N^~ ot {q = 0) is larger 
than that of the smooth signal and consequently, this quantity is a criterion for the roughness of 



9 



15 I — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




250 500 750 1000 

i 




a 



FIG. 4: Upper panel corresponds to generated series for various values of roughness exponent (from top to 
bottom: H = 0.8, H = 0.6, H = 0.5 and H = 0.3). Lower panel shows the level crossing analysis of time 
series with above mentioned Hurst exponent. Results regarding the original and shuffled series have been 
indicated in each panel. 



the processes. For a correlated signal, N^ t (q = 0) is smaller than that of a completely random one, 
while for anti-correlated data set N^ t (q = 0) is larger than the same series for which its memories 
have been destroyed by shuffling. Figure (j3J) shows the level crossing results for some fractional 
Gaussian noises (fGn) and the corresponding shuffled signals. For an anti-correlated series the 
so-called Hurst exponent is H < 0.5, for a completely random Gaussian data set H = 0.5 and for a 
correlated signal H > 0.5 j^. In addition, for moments q < 1, those terms in N^ t (q < 1) will 
become dominant that have small deviations from mean level, a, which demonstrate the statistics of 
small fluctuations. On the contrary, for q > 1, those terms with large fluctuations in the integrand 
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FIG. 5: Algorithm of map simulation and analysis used in this paper. 



of Eq. (|lip are dominant in N^ t (q > 1) explaining the statistical properties of the upcrossings far 
from mean level. These terms correspond to the rare events. 

In what follows, we are going to derive v£ for discrete temperature fluctuations in the simulated 
map with and without cosmic strings on the flat sky and compare them to find a robust criterion 
to distinguish pure Gaussian fluctuations from fluctuations containing cosmic string components. 
The motivation of using this simple method can be justified according to the following reasons: 



i) Since the statistical isotropy is valid as a major statistical property 6814701]. one can cut many 
one-dimensional signals in every directions and by ensemble averaging, compute the frequency of 
upcrossings (downcrossings) at all interested levels along with their variances. 

ii) As one can see in the next section, the superposition of fluctuations produced by cosmic strings 
can generate new and extra ups and downs in temperature fluctuations, consequently finding such 
statistically meaningful footprints in the map in comparison with pure Gaussian signature including 
instrumental noise may potentially help us to get deep insight in the cosmic string detections. 

iii) Due to the phase coefficient in the Fourier analysis, it seems that many trivial imprints of 
cosmic strings diminish or at least are mixed with other observational phenomena so it is another 
motivation to investigate the imprint of cosmic strings in the real space as mentioned in Ref. ^ 



III. SIMULATION OF MOCK CMB MAP 



In this section, we describe in detail our sate-of-the-art code which has been written in Fortran 
an guag e by authors to simulate temperature fluctuations maps (also for more details see refs 
19-48,0, Q)- Generally, program has four parts: The first part creates pure Gaussian 
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FIG. 6: A Gaussian map for ACDM model based on WMAP-7 observations. Map size is 10° with R = 1.5'. 



fluctuations corresponding to the standard inflationary model in the presence of various models 
governing the evolution of background. In this paper we use ACDM model in the flat Universe. 
Nevertheless, our program can be easily modified to other cosmological models governing the 
background evolution of universe such as quintessence model and so on {73! ]. The second part 
contains anisotropies produced by long cosmic strings by means of Kaiser-Stebbins effect. The 
superposition of Gaussian and cosmic strings anisotropies as well as the expected instrumental noise 
are produced in the third part. The normalization factor for each component will be computed in 
this part. Finally, the algorithm for searching the signatures of cosmic string relying on the level 
crossing analysis is given in the last part. To find more reliable results and reduced statistical 
errors, an ensemble averaging over at least 200 runs for a set of input variables will be done. As 
one can see in the following, we need a fine resolution map. However the random behavior of 
cosmic strings limits our simulation on the flat sky map approximation In this case 

the Fourier modes are bases functions instead of being spherical harmonics. Before going further, 
for pedagogical purposes, we summarize all procedures to be done for map making in Figure ([5]). 



A. Gaussian map on the Flat sky 

In this subsection a Gaussian map for temperature fluctuations will be generated. We use 
ACDM model and the best fit values have been inferred based on the most recent observations 
such as WMAP-7, Supernova type la gold sample (SNIa) and Sloan Digital Sky Survey (SPSS) 
with the most familiar initial power spectrum established by standard inflationary scenario [76. \t\. 



12 

As discussed in details in Refs. Q^Q] and indicated in Figure ©, first of all, we should determine 
the values of initial parameters relevant to the Gaussian map. A part of these parameter are related 
to the cosmological frame-work and others are the size of simulated map, G, angular resolution, 
R, and finally the variance of instrumental noise, namely cr n oise- Since we are interested in flat 
sky, namely a simulated map with size less than 60°, two-dimensional plane wave (Fourier basis 
functions) are used instead of spherical harmonics. The stochastic temperature fluctuations field 
on the flat sky will be written as: 



^r(k X ,ky) = ^^^l[Z 1 (k x ,k y )+iZ 2 (k x ,k y )} (12) 

Actually to guaranty the Gaussianity of distribution based on central limit theorem, we add 
Zi(k x , k y ) and Z2(k x , k y ) as two independent Gaussian random numbers with zero mean and unit 
variance. The factor 2 in the denominator stands for canceling the corresponding variance of Gaus- 
sian random numbers. To set the value of Ci for a given set of (k x , k y ) in inverse degrees, one can 
compute I as: / = ^■■\/k% + k^. The initial power spectrum, C;, is determined by running CAMB 



software 



46|, 



781 ] . The power spectrum for non-integer values of / can be evaluated by linear inter- 



polation in the produced C\. To this end we modify the public CAMB software for underlaying 
theoretical model and set the best fit values of corresponding cosmological parameters based on 
constraints given by the most recent observations. Finally to keep the statistical isotropy and to 
diminish the undesired preferred direction appearing in simulation, we construct four independent 
Gaussian maps and superimpose these separate components according to 

AT 



47|]: 



— {x,y) G = 

1 [AT , AT . AT , AT 

2 -jr{X,y)Gl + y^mH - X,y) G 2 + — ~ V)G3 + yl^max ~ ^, 2/max ~ VjGi 



13) 



The pre-factor 1/2 is necessary to keep the initial standard deviation. As mentioned before, in 
this paper we simulate Gaussian temperature fluctuations based on ACDM H] including best 
fit values for cosmological parameters based on WMAP-7 observations. Figure ([6]) shows a typical 
simulated Gaussian map for ACDM model. 



B. Fluctuations Produced by Cosmic Strings 



In this subsection we explain in detail the algorithm used for generating cosmic strings by using 
the Kaiser-Stebbins effect which originally discussed by Perivolaropoulos 38j ( for more details 
see Refs. 38, 46-48, 0, Q]). In our toy model, since the size of a typical loop of string is about 
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z+i 



FIG. 7: Typical size of loop in arc-second units as a function of redshift. 



79 



80( | (see Figure ([7])) an d the resolution of simulated map is of the order of one 



arc-minute so we can ignore the contribution of these closed strings and suppose that the main 
anisotropy in the string map is due to the straight cosmic strings. In addition, it has been shown 
that loops contribute to temperature fluctuations via the well-known Sachs- Wolfe effect and due 
to the small bandwidth of simulated map the overall effects can be represented as Gaussian noise 
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8l|. 



The scaling behavior of the correlation length scale of straight strings demonstrates that the 
number of strings crossing a given Hubble volume to be fixed. In another word the number density 

Hit 
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80 



86]. Therefore, one 



of strings to be invariant under rescaling of correlation length scale 
can deduce the numbers of straight strings crossing each Hubble volume at any given time to be 
M = 10. Their overall properties such as orientations and velocities are statistically uncorrelated 
for length scale larger than the corresponding horizon scale 



79 



80, 



87j. The time interval in our 



simulation is limited to the interval [t/ 5 ,to] an d since the cosmic strings have relativistic velocity, 
one can suppose that after 2tn {pH is Hubble time), a new Hubble volume will be generated and 
consequently a new network of cosmic strings affects the propagation of the photons transmitted 
freely after last scattering surface based on the Keiser-Stebbins phenomenon. This is also justified 
regarding the fact that a photon ray travels most of the Hubble volume during twice the Hubble 
time, and after that this ray encounters a new volume and will be affected by a new string network. 
So one can suppose that at successive time intervals according to the relation, tj+i = 2ij, the photon 
ray could be affected by a new string network. According to this statement, the number of separate 
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FIG. 8: Left panel shows a sketch of an extended window used to simulate the effect of one cosmic string 
on the CMB temperature of 6 x 9 map at redshift z,-. Generally, strings located in the extended map of 
size (G + 20,) 2 can affect the simulated patch. Right panel corresponds to a schematic view of horizons 
at successive steps to accumulate KS effects on the temperature fluctuations through traveling from the 
last scattering surface toward an observer. After about two times of the horizon length scale, photons are 
encountered with a new Hubble volume, consequently, new string network makes kicks on them. 



string networks (patches) to be simulated is given by (see Figure ([8])): 

to 



N P = log 2 



15 



(14) 



In the string map for the recursion condition, = 2ij, we let the new network of cosmic strings 
to kick the photon paths. The apparent angular size of horizon at Zi is given by: 

Jz, 



e, 



H(z;{H}) 



-T 



VI SZ ^Uo H(z;{£}) 



(15) 



Hq^/Wk\ 

here J-(x) is x, sin(a;) and sinh(x) for flat, closed and open universe, respectively. In the above 
equation, H(z; {£}) is the Hubble parameter and {£} includes all model parameters. The number 
of straight cosmic strings with size 6 t - at Zj which should be simulated and placed randomly on 



the desired simulated map will be 



46 48, 



2^: 



Pi = M 



[6 + 20j 



(16) 



where is the size of simulated map with resolution equal to R, consequently the number of pixels 
for extended map at ith step of simulation is [0 + 2@i]/R. Actually, to retain all fluctuations 
produced by cosmic string in a limited area of sky which has been simulated we establish an 
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z+1 z+1 

FIG. 9: Left panel corresponds to the number of straight cosmic strings in unit of M (M is the number of 
cosmic strings at each Hubble volume given by the scaling solution) as a function of redshift for a simulated 
map with size of 2.5° x 2.5°. Right panel indicates cumulative string number as a function of redshift for 
various values of simulated map size for a flat Universe with O m = 0.27 and Q\ = 0.73. 



extended area for each step. The size of mentioned extended window will be + 2@j at ith patch 
and our simulated map to be considered at the center of this window. Left panel of figure (|8]) 
shows this configuration for fluctuations produced by Kaiser-Stebbins phenomena on the cosmic 
background radiation fluctuations. Upon the calculation of the size of horizon for ith successive 
step by using Eq. (fT5j) . we embed a straight string and apply a shift for temperature fluctuations 
according to Eq.([T]). Since the mean value of fluctuations in the simulated map has no physical 
significance, we add half of (AT/T) to one side of string and subtract the remaining value from 
other side, namely we use 

(J^r) =±AvGn ls v s r (17) 

where r = |n.(v s x e s )|. The direction of observer for a small field in the map is approximately 
constant. Since the string's velocities and orientations are random and furthermore r takes into 
account the contribution of orientation and projection effect of cosmic string, consequently its 
value will be given by generating a uniform random number over the interval [0, 1]. In addition a 
binary flag is produced to decide which side of string to be positive and which one will be negative 
for temperature fluctuations. The location of string in the extended window is also demonstrated 
by an additional proper random number. We take into account the contribution of projection 
from 3-dimensions to 2-dimensions of cosmic string by multiplying the proper size of strings at 
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FIG. 10: A string map for Gfi = 1 x 10~ 7 . Map size equates to 10° with R = 1.5'. 



each patch by the cosine of the angle deriving from [0, ^] interval randomly. We fix the relevant 
coefficient for string simulation e.g. the number of cosmic string for each Hubble volume, M = 10 
and Vg'jg = 0.15 [88l |. It must be pointed out that the fixed values for string quantities change 
significantly from one theoretical model to the other, consequently in this toy model, we take most 
popular values for mentioned relevant quantities {J. As shown in Figure Q the number of strings 
for each step depends on the size of simulated map as well as the underlying cosmological theory. 
Finally, temperature fluctuations simulated for various patches should be superimposed to reach 
the final situations for a given pixel in the simulated map. Figure ()10p displays a string map for 
which the cosmic tension is G\i = 1 x 10~ 7 . As an additional check for consistency of our simulation 
we compute the power spectrum of our pure string simulated map for different values of G\x as well 
as pure Gaussian generated map in small angle approximation 0, [|l|, Q] . Figure (fTT|) indicates 
the angular power spectrum for pure Gaussian and pure cosmic string simulated maps for typical 
field of view of 2.5° and a resolution of 1.5'. Since the size and resolution of generated map is 
finite and temperature fluctuations are limited to the flat sky, the multiple moment will run from 
I ~ 10 2 to 10 4 . In addition the expected scaling behavior for pure cosmic string map satisfies the 
relation 1(1 + l)Ci oc l~ v for I S> 1 with rj ~ 0.901q'q5. Let us explain the final task regarding 
the simulation of a Gaussian map including induced cosmic string and instrumental noise. As 
discussed before we simulate all components separately and according to the linear perturbation 
theory, we superpose all sources to include all phenomena and reach the final situation. We do 
not expect the observed temperature power spectrum to be changed in the presence of cosmic 
string, consequently the multiplication by a factor, u less than unity causes the power spectrum to 
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FIG. 11: Angular TT power spectrum of CMB fluctuation. Solid line is calculated by CAMB software 
for the best fit values of ACDM based on 7-year WMAP data. Dashed line corresponds to a simulated 
Gaussian map with size 2.5° x 2.5° and map resolution is R — 1'. Long-dashed, dashdd and dashdot lines 
correspond to fluctuations Long-dashed, dashdd and dashdot lines correspond to fluctuations generated by 
cosmic strings for some typical values of G/j, indicated in the plot. The power-law fitted function is also 
presented in figure with the scaling exponent equates to r\ — 0.90^o'gg for I ^> 1 [3l[. 



remain unchanged. The instrumental noise does not modify the real sky map and its significance is 
demonstrated by the maximum amplitude of noise represented by 1 . This maximum will be 

J 1 1 J 1 \max 

taken as = WfiK throughout this paper according to the anticipated instrumental noise 

in South Pole Telescope (SPT) 891 ] . To finalize our simulation, we use the same method explained 



in Ref. [461 ] . The power spectrum of simulated map can be read as follows: 

Cl(G+S) = ^ 2 Q(G) + Cl(S) ( 18 ) 

The cross-correlated term is zero due to independency of Gaussian and string components. We 
demand the left hand side of the above equation to be equal to power spectrum derived from 
observation, given by CAMB software. Therefore for each / from l m i n to l max , we find a separate 
cj's. To compute a single cj, we rely on Bayesian statistics for observations, {X} : {w(/)} and model 
parameters, {Q}. We define Q: X 2 (ty = Tfi=i max [w o bs(0 — w theG; ^)] 2 - Here w o b s (0 is computed 
directly from Eq. (jl8p and wthe(^;^) derived by an arbitrary merit function which is supposed to 
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FIG. 12: Left panel shows the Gaussian map. The middle panel corresponds to pure cosmic string for 
G/i = 1 x 10 . The superposition of all components is shown in right panel. The resolution of these map 
is R = V. 
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FIG. 13: Left panel shows the Gaussian map. The middle one corresponds to Gaussian with cosmic strings 
added for G/i = 1 x 10~ 7 . The superposition of all components with instrumental noise is shown in the right 
panel. Maximum value for temperature fluctuations created by noise is considered to be (^f~~ max ) N = 10(iK. 
The size of all maps is 2.5° with a resolution of R = V . 



be a first order polynomial function throughout this paper. The y- intercept of merit function will 
be assumed as a best fit value of oj [46]. It must be pointed out that to improve the reliability 
of u), we do various simulations and finally by ensemble averaging, an averaged value for uj can 
be retrieved. Figure (|12|) indicates different components simulated without expected instrumental 
noise. Figure (|13p shows simulated map including instrumental noise. Now every things ready to 
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FIG. 14: Left panel shows the 1 + 1-dimensional increment of temperature fluctuations for Gaussian and 
superimposed cosmic string component with G/i = 2 x 10~ 7 . Right panel corresponds to the same as left 
panel except including instrumental noise. 

investigate temperature fluctuations on the fiat sky with and without cosmic strings. We rely on 
a robust statistical method namely, level crossing analysis to explore the capability of our method 
to detect the cosmic strings for various values of string's tensions. 

IV. ANALYSIS OF CMB MAP 



As explained in details in section II, we are interested in upcrossings the temperature fluctuations 
at an arbitrary level, a. Or according to probability point of view, we should compute the volume of 
joint probability density function which is satisfied in the proper conditions mentioned in section 
II (see Figure ([3])). We expect the fluctuations in the presence of cosmic strings to be rougher 
than pure Gaussian fluctuations. Since the level crossing analysis is applicable for one-dimensional 
series we construct 1 + 1-dimensional series from simulated map. Furthermore due to the statistical 
isotropy property of simulated map 



68H70I]. the 1 + 1-dimensional data sets will be created in all 
given directions as different ensembles for analysis. Actually these series contain temperature 
fluctuations as a function of pixels in one dimension in an arbitrary direction. Hereafter we turn 
to investigate 1 + 1-dimensional signals as an input for level crossing subroutine which is the final 
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part of the main simulation program. We introduce the increment of fluctuation as: 

AT(k) s _ c r '*-')-n (i9, 



T 

In Figure (|14[) we have plotted a typical 1 + 1-dimensional increment series for Gaussian and 
cosmic string signals embedded in Gaussian fluctuations. We see that in the presence of strings, 
fluctuations near the mean as well as far from the mean have been changed. The magnitude of 
these changes depends on value of cosmic string tension. As pointed out before, we enumerate 
the number of crossings with positive slope at whole existence level for all 1 + 1-dimensional data 
set in a map for a given G\x. It is worth noting that according to conservation law, for long-run 
simulation all crossings with negative slope are statistically equivalent to that of with positive slope 



60J. Based on this analysis, we find many quantities which may be potentially used as a criterion 
to distinguish between the maps with and without straight cosmic strings. In Figure f)15|) we plot 
the u£ versus a for various values of Gu. Generally, total crossings in the Gaussian map is less 
than that of with cosmic string. As we expect, the existence of instrumental noise decreases the 
efficiency of our analysis. To compare simulated maps we refer to the generalized form of roughness 
of signals which is given by Eq. . N^ t (q). The various values for g's take into account crossings 
with different measures in the total level crossings. Figure (I16p indicates the generalized roughness 
function for some typical values of Gfi. To compare the value of generalized roughness function we 
apply Student's t-test for equal sample size and unequal mean and standard deviation for each q's 



which is defined by: 



t(q) = (N+ t (®,q)-N+ t (®,q)) x ^ 



2 ,f r '"" : (20) 

here the symbols © and ® stand for Gaussian (G) and Gaussian + String (GS) for simulated map 
without noise, respectively. For simulation including instrumental noise these symbols are replaced 



by Gaussian + Noise (GN) and Gaussian + String + Noise (GSN), respectively. The N^ t is 
the ensemble average of generalized roughness function for each q's. The corresponding standard 
deviation is given by cr(q). N run indicates the number of simulated ensembles. The P-value, p(q) 
for t(q) based on t-distribution function with 2N run — 2 degrees of freedom, will be determined. 
By the combination of above P-values we introduce the following quantity: 

Qmax 

X 2 = ~2 Y, ln ^) ( 21 ) 

Q — Qmin 

Finally, the last P-value, Pfi na i associating to the x 2 by using chi-square distribution function 
for 2 ( i max ~Q min j _ 2 degrees of freedom will be computed. For 3a significance level namely 
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FIG. 15: Level crossing for increment of temperature fluctuations for various simulated components. The 
value of G[i for left and right panel are 1 x 10~ 6 and 1 x 10 -7 , respectively. 



P final < 0.0027, we can conservatively say that there exists a significance difference. According to 
the value of P final indicated in Figure (|17p . we conclude that our method can detect the cosmic 
string with Gfi >4x 10~ 9 and in the presence of anticipated instrumental noise the lower bound 
is G\i > 5.8 x 10~ 9 . The final strategy for detecting cosmic string is as follows: At first by using a 
map from observation we compute and N^ t (q). Second we simulate a pure Gaussian map with 
expected instrumental noise and again the level crossing and generalized roughness functions to be 
calculated. The P-value for deciding about the existence of significant difference between N^ t (q) 
for real and simulated Gaussian maps will be determined. In the case of significant difference, 
we change the value of Gfi for simulated Gaussian map including cosmic strings and compute the 
corresponding P final- Upon the P final for simulation to be equal to that of given for observation 
we conservatively can express the existence of cosmic strings in the observed map with mentioned 
Gfi. 

V. SUMMARY AND CONCLUSIONS 

In this paper we relied on a robust method in complex system which has been introduced 
originally to investigate the statistical properties of 2 + 1-dimensional surface by mapping to 1 + 1- 
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FIG. 16: Generalized roughness function versus moment, g, for Gaussian+ String + Noise and Gaussian + 
Noise. The size of symbols are almost equal to their error bars. Here we run over 200 ensembles. The value 
of string's tension from left to right are G/j = 1 x 10~ 7 , G/i = 1 x 10~ 8 and Gfi = 6 x 10 -9 , respectively. 




FIG. 17: P-value as a function of G/i. We run over 200 ensembles to compute each point of this plot. To 
make more obvious we fitted data by typical fitting functions. 
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dimensional data sets 60-65]. Using this method we explore the capability of finding the footprints 
of straight cosmic string based on the KS effect on the cosmic microwave background fluctuations 
since the last scattering epoch up to now. To this end we take a simple toy model to simulate the 
effect of cosmic strings demonstrated in Refs. j^E^]. The contribution of loops in this model can 
be ignored (see Figure (0). The superposition of kicks due to the embedded cosmic strings cause 
some additional statistical fluctuations in simulated map. Consequently there exists a difference 
between the number of level crossing for whole reachable level in f + 1-dimensional fluctuations of 
pure Gaussian map in comparison with map containing cosmic strings. To quantify this deviation 
we introduced an increment signal and use robust quantity to take into account all deviations 
relative to pure Gaussian fluctuations, namely generalized roughness function. For q = the 
generalized roughness function corresponds to sum of upcrossings at all fluctuation levels which is 
a measure of roughness of underlying signals. Based on joint probability density function of signal, 
p(y = a,y'), we found that the fluctuations near the mean fluctuations could be a proper criterion 
to discriminate pure Gaussian map from the map induced by cosmic strings. Our result shows 
that level crossing analysis could place a bound for cosmic string tension G[i >4x 10 -9 for 200 
simulated maps with R = V . To explore the contribution of instrumental noise we relied on South 
Pole Telescope experiment and took (~^\ max ) N = WfiK for noise map. The new bound for map 
including instrumental noise is G[i > 5.8 x 10 -9 . The small change in the bound on Gfj, due to the 
noise could be justified regarding Figure (|15[) . The sensitive part of v a is crossing near the mean 
value so the presence the undesired noise can only have minor statistical effects on the ability of this 
model to place a lower bound on the string's tension for small value of G/x. In addition, our results 
for both simulations were not affected by increasing the size of ensemble up to 200 corresponding 
to increase the coverage of sky map. For the realistic observations there are additional problems 
to be addressed. The first one is the foreground and the instrumental noise which is supposed 
to be considered by adding expected noise and smoothing our simulation. As mentioned in Refs. 



46H48I] discontinuity related to the scanning strategy also puts extra imprint on data set. It must 
point out that our method is sensitive to overall fluctuations in all directions, therefore, since it is 
expected that by increasing the coverage of underlying map the statistics of fluctuations near the 
mean to be significant, the mentioned spurious phenomenon will have a negligible contribution to 
the final result. However we can remove all parallel discontinuities in preferred direction imposed 



by scanning maneuvering of observational instrument 



use more realistic models 
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84, 



46]. Final remark is that we would like to 



911 ] to simulate a map taking all contributions of cosmic strings on 



the cosmic microwave background into account and implement our method to explore the effect of 
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cosmic strings in our future works. In addition as a future study, it is useful to apply this method 
to distinguish the effect of generic cosmic strings and superstrings on the temperature fluctuations. 
The n-point peak-peak correlation function is also another method which could be used for the 
same purpose. 
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